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THÈSE  D’A^ALYS 


RECHERCHES  SUR  QUELQUES  FORMULES  D’ANALYSE, 
PARTICULIER  SUR  LES  FORMULES  D’EULER  ET  DE  STIRLING 


Ce  travail  est  divisé  en  deux  parties.  Dhns  la  première  partie,  je  fais  con- 
naître Forigine  et  la  loi  des  nombres  de  Bernoulli  ; dans  la  seconde,  j’établis 
d’une  manière  rigoureuse  la  formule  d’Euler;  je  donne  l’expression  du  terme 

complementaire;  j’en  déduis  ensuite  comme  conséquence  la  formule  de 
Stirling. 


PREMIÈRE  PARTIE. 

I. 

ORIGINE  DES  NOMBRES  DE  BERNOULLI;  EXPRESSIONS  GÉNÉRALES  ET  PROPRIÉTÉS 

DIVERSES  DE  CES  NOMBRES. 

1.  Les  nombres  remarquables  dont  nous  nous  occuperons  dans  la  première 
partie  de  ce  travail  sont  connus  des  géomètres  sous  le  nom  de  nombres  de 
Bernoulli  ou  nombres  Bernoulliens.  C’est  effectivement  Jacques  Bernoulli  qui 
les  a rencontrés  le  premier  en  étudiant  les  séries  dans  lesquelles  se  développent 
certaines  fonctions  simples,  et  depuis  ils  se  sont  présentés  dans  un  grand 
nombre  de  questions  diverses. 

Pour  définir  ces  nombres,  considérons  la  fonction 


e‘—i 

dans  laquelle  e désigne  la  base  des  logarithmes  népériens. 


r 


Cette  fonction  se  réduit  à l’unité  quand  on  fait  z = o ; on  peut  donc  écrire 


= 1 4-  £, 


c' — I 


désignant  une  quantité  qui  s’annule  avec  z;  on  tire  de  là 


d’où 


I z -h  I — e* 


S 

z 


I 

I .2 
I 


Z 


l , 2 


On  voit  que  pour  z = o,  ^ se  réduit  à — et  si  l’on  pose 

E 1 , 

- = h cpz, 

z 2 ‘ 


la  fonction  9 z s’annulera  en  même  temps  que  z;  on  a 


La  fonction  <pz  est  impaire;  en  d’autres  termes,  elle  change  de  signe  quand 
on  change  z en  ~ z.  En  effet,  on  a 


donc 

Z)=z-(pz. 

La  fonction  <pz  reste  finie  pour  toutes  les  valeurs  de  z réelles  ou  imaginaires 
dont  le  module  est  inférieur  à 2 7:  : donc,  d’après  un  théorème  célèbre  de 
M.  Cauchy  (* *),  cette  fonction  sera  développable  en  une  série  convergente 


(*)  Ce  théorème  a'été  énoncé  pour  la  première  fois  dans  un  Mémoire  publie  a funn  en 
i832  par  M.  Cauchy.  On  lit  effectivement  dans  ce  Mémoire  : 

« La  fonction  f {x)  sera  développable  en  une  série  convergente  ordonnée  suivant  les  puissances 

* de  X,  si  le  module  de  la  variable  réelle  ou  imaginaire  x conserve  une  valeur  inférieure  a 
» celles  pour  lesquelles  la  fonction/ [x)  cesse  d’être  finie  et  continue.  » 


( 5 ) 

ordonnée  suivant  les  puissances  ascendantes  de  z,  tant  que  le  module  de 
cette  variable  sera  inférieur  à an.  Et,  parce  que  la  fonction  oz  est  impaire 
on  aura  un  développement  de  la  forme 


— -B.2 


B, 


1.2  2.3.4^^*,.2.3.4.5.6--‘ 

Les  coefficients  B„  B„  B3, . . . , B„,  sont  dits  les  nombres  de  Bernoulli. 


Développement  de  cotx  en  série  ordonnée  suivant  les 

croissantes  de  x. 


puissances 


2.  Pour  obtenir  ce  développement,  nous  partirons  de  la  formule 


I I 1 


2 e*— I Z 

que  1 on  peut  écrire  de  la  manière  suivante  ; 


I e'i  -\-e  2 , 

QZ  ~ 

^ 2 £ _£  Z 

ei  — e 2 


en  posant 
il  vient 

d’où  l’on  tire 


Z=z  u\J  —i, 


1 «2 

2 


-y/—  -tJ~ 

•1'  I „ r,  V ‘ 


e 2 


cot  — U 
I 2 


V'-'  usj~) 


V — > U \J~i 


U ~ 2«^ot2“--V-i-9(«v'-i)=B,  — + B, — 


*•2^  “^i.2.3.4  1.2. 3.4.5,  e'*""’’ 

formule  qui  subsiste  pour  toutes  les  valeurs  réelles  ou  imaginaires  de  u dont 
le  module  est  inférieur  à 2 tt.  Faisant 


U — 2 JT, 


il  vient 


( 6 ) 

I „ 2’.Bj  3 . 2’.B<  - 

(i)  cotx— aBjjr  'ÿj  ^ ~ 34.5.6^  34.5.6.7 

mais  cette  nouvelle  formule  ne  subsiste  que  pour  les  valeurs  de  x dont  le 
module  est  inférieur  à n. 


Développement  de  tanga:;. 


5.  Si,  dans  la  formule 


Z Z Z B|Z’  BjZ* 

— --  = ! h Z(Ù  Z = \ 1 5-7 

— I 2 ^ 2 1.2  1.2. 3. 4 


nous  changeons  z en  2 z,  il  vient 


2 Z 2 Z 2’.  Bi  z’  2*. B,  Z* 

— — — I + 5 — T + . . . . 

e-‘  — I 2 1.2  1 . 2 . 3 . 4 


On  a d’ailleurs  identiquement 


+ 1 e- — I — i ' 


donc 


Z Z B,  z‘‘  B’ Z*  r 2z 

e‘4-t  2 1.2  1.2. 3. 4 L ^ 


OU 


OU  encore 


2’B,  z’  2‘B,Z‘  -f- 

1.2  12.3.4 


B2Z' 


1 I /O  M3  2 / * \ 

— — = (2-  — i)B, (2*  — I) 5-7  + ...; 

2 e--\-I  ’ I .2  ^ 1.2. 3. 4 

changeant  z en  2 z et  remarquant  que  l’on  a identiquement 


1 c*  — I 1 — e 2 


2 e’  -f- 1 2 e*  + 1 2 


e2-+-(?  2 


1 c*  — e~‘ 

2 + 


= 2 .(2“  — l) 


B^ 

1.2 


2*  (2*  -1)63-- 


.2.3  4 


il  vient 


Si  l’on  fait  maintenant 


( 7 ) 


Z = x\j~l^ 

on  pourra  écrire 

(2)  tanga:=2*(a»-i)B,^-2*(2*-i)B»-^! H 

' 1.2. 3. 4 

formule  qui  ne  subsiste  que  pour  les  valeurs  de 

rieur  à -• 

2 


2«(2“-i)B3 


2.34.5.6  •” 
X dont  le  module  est  infe 


Développement  de  coséco:. 


4.  Nous  avons 


coséc  X — =:  ^ ^ 1 rnt  ' .r-  _i_  ' t ' 

sm.r  2sin|xcosix  ~ 2 2 ~ 

^coti2r=i-B,f_B,^-B 


2 

I . I 


'2  "2.3.4  2.3.4.S.6  •••’ 

ï34 


jtang,.  = - ,)  B,  , + (a-  - ,)  + (a-  - 0 B. 


donc 


(3)  coséc  X — - 0 + (2® l) 

^•4  3. 4. 5. 6 


II. 

LOI  DES  NOMBRES  BERNOULLIENS  ; EXPRESSIONS  DIVERSES  DE  B,.  EN  FONCTION  DE  «. 

O.  L’expression  générale  des  nombres  de  Bernoulli  a été  donnée  par  La- 
place  ( rojez  le  Traité  de  Calcul  différentiel  et  intégral  de  Lacroix 
tome  III,  page  107).  Mais  la  méthode  suivie  par  l’illustre  géomètre  est  très- 
compliquée,  et  nous  croyons  pouvoir  présenter  ici  une  solution  plus  simple 
et  plus  élégante  de  la  question. 

On  a 

t sin  .r  . , i 

tang  X — = sin  x (i  — sin^  x'f  2 , 


( 8 ) 

_ j_ 

et,  en  développant  (i  — sin^x)’  2 par  la  formule  du  binôme,  il  vient 


, . I . 3 3.5  . 3.5.7...  2 « — 3 . 

1 tanc  JC  sin  .7:  — sin^  jcH 7Sin^  jc  h 7— ^sin’  x + ...  -t sin"‘“‘ jc 

\ . ^.4  .,4,e 

I ^ 3. 5. 7-. .(2/2 — 3)f2« — 1)| 


2.4.6...  (2  « — 2).2  rt 


2.4.6  ’ — 2 

[sin^"'*"'  JC  + etc.  ]. 


Le  coefficient  de  æ‘"~'  dans  le  développement  de  tang  x en  série  ordon- 
née suivant  les  puissances  de  x est 


î (‘Il 

1.2...  (2«  — l)  \ 


‘ tangx 


dx‘'‘~  ' 


d’après  la  formule  de  Maclaurin  ; on  a donc 


'^’""'-tangx\  _ (2«  — 1)(2«  — 2)...3.2^2«/„2«  .^T3  — Op 

, jo~~  2.3.4...2W 


OU,  à cause  de  la  formule  (4), 


[ Id^’'  ' tang.r\  sinjc\  i (t 

V dx^'>-'  V dx^'^-'  i V 


(5) 


d^"~'  sin’  X 


3 /c?*"”' sin^jc 


dx™~'  )q  2.4  \ dx-''  jQ 

3.5.7. . .(2rt  — ^ 

2.4.6.. .(2rt  — 2)  \ C?JC^"~ 


mais  le  développement  peut  être  arreté  au  dernier  terme  écrit,  car  il  est  évi- 
dent que  pour  .r  ==  o,  la  dérivée  d’ordre  2 « — i de  la  fonction  sin  j:*"-^-*-** 
se  réduit  elle-même  à zéro.  Cette  fonction  est  effectivement  développable  en 
série  ordonnée  suivant  les  puissances  ascendantes  de  x,  et  le  degré  de  son 
premier  terme  est  2/2  — i 4-2^:;  en  sorte  que  la  dérivée  d’ordre  2/2  — i con- 
tiendra le  facteur  x dans  tous  ses  termes. 

Gela  posé,  différentions  2 22  — i fois  la  formule  connue 


sin  x= 


J2W— l)(2W-2)...(/72-t-l) 

1.2.3.  ../K  1 1.2.3.  ..(w  — 2) 

(2m  — i)  (2.  m — 2).  . .(m  — 1)  . ►. 

- ‘ sin  5 — . . . 

1 . 2 . 3 . . . w — 3 

{2m — 1)  ( 2/w  — 2).  . .f/w — f'4-2)  - / . \ 

— — 5 ^ sin (21  — ijxzfz... 

I .2.3 . . .{m  — i)  ' ' 

± sin  (22/2 — i).r. 


nous  aurons 


( 9 ) 


(2m  — i){2m  — 2).  . .(m  -|-i) 

rTTf^rrT] ~ COS  JC 

__  ^2«-l  — 0(2W  — 2).  . ./w 

1.2.37  .(/«-  2) cos3jc 

(2W  — l)(2w  — 2).  ■ ,f/«  — 

1.2.3.  . .{m  — 3 j COS  5 JC 

:h(2/n  -i)'*”-*cos(2/72-i)x;  ‘ 


posant  successivement 


m 


et  faisant  en  même  temps  jc  = o,  il  vient 

f s'in.7;\  , 

\ ; 0 (~  0 

-1 

V /o  2»  Ll  ~ ^ J’ 

_ ( — 0"'^'  ['5.4  5 

2*  L I . 2 I 


( 


sin’ . 
dx‘"~  ‘ 


f\  [1-6. 


[tTI  - ril  3="-' + i 5“-' - , 


sin2«-i  ^ — I }"+' 

. dx‘‘'‘~' 


lif j)  (2«  — 2j.  . (2«  — l)(2«  — 2).  . .(/z-4-2) 

1 .2.3.  . .A<  — J 1.2.3.  ^ O " * 


. n 


{%n  — i)  . 

1_  — I — -(^2  n — 3 

±:  (2  n — ] 


11  nous  est  aisé  maintenant,  au  moyen  de  la  formule  (5),  d'écrn-c  l'expres- 
sion  e ^ ^j.2n-i  6t  comme  cette  expression  est  égale  à B„.^ — 


2 


on  pourra  poser 


( lo  ) 


ryin  / o2« 


2’"  (2^"  — i) 


J 

.^2Li  J 

1.3  p5.4  ^‘^2/1-1  , 1 

i^2.4Ll.2  I J 


1.3.5  r 7.6.5  7 32«— < _|_  2 52"-<  «2.-!-) 

Li.2.3  1.2  I '• 


2®.  2.4.6 


I .3 .5.  . . (2«  — 3) 
2^"—^.  2.4.6.  . . (2  rt  — 2) 


(2/2  — l)(2/î- 

- 2). ..(«  + .) 

1 

1.2.3.. 

.(»  — i) 

1 

( 2 « — 1 ) ( 2 « ~ 

-2}...n  + 2 

1.2.3.. 

.«-2 

1 ^ 

(2«  — l)  - 

^ ( 2 n 

- 3)^"-' 

U 

(U  fl 

? 

telle  est  l’expression  des  nombres'^Bernoulliens  que  Laplace  a obtenue  le  pre- 
mier, en  suivant  une  marche  très-différente. 

La  manière  d’employer  cette  formule  est  évidente.  Pour  trouver  B,,  on  fera 
n = J , en  s’arrêtant  à la  première  ligne  ; pour  trouver  on  s’arrêtera  à la 
deuxième,  en  faisant  n = etc. 

On  obtient  ainsi 


83  = 7^, 

* 42 


B. 


66 


- (n- 


6.  Avant  d’aller  plus  loin  dans  l’étude  des  nombres  Bernoulliens,  nous 
démontrerons  une  formule  dont  nous  ferons  un  fréquent  usage. 

Cette  formule  est  la  suivante  : 


(6) 


I I I 

cot-  — 

VIT. 


dx. 


(’“)  M.  Schlomilch  a publié  sur  ce  sujet  un  travail  que  j’ai  consulté  utilement. 


( ) 

Pour  I établir,  partons  de  la  formule  connue 


ï 1 


COtx=:)^  ^ ^ — 2 7T  X — 3it 


+ ... 


' H 


jt  + tt  X ~ 2n  x-i-3, 


Cette  série  est 
— 7T  et  -f-  ;r. 


convergente  pour  toutes  les  valeurs  de 


X comprises  entre 


En  effet,  on  peut  écrire 


cotj:  — - 

X 


I 

2t:  ■+■  X 


nn  + x 


» 


OU 


COt  JT =: 

X 


— '2  JC  [ î 

L’ï'*  — x^ 


4 7r’  — X- 


9 7t’ X^ 


«’tt* — x^  • . .J- 


r e développcmern  étant  mis  sons  cette  forme,  on  vo.t  que  pour  des  valeurs 
de  ,z  comprises  entre  - tr  et  + n,  les  termes  entre  crochets  à partir  du 
deuxieme,  forment  une  somme  moindre  que 


pan 


[3  8 75 

Or  cette  somme  a une  limite. 

Cela  posé,  considérons  l'intégrale  définie  = i,  qui  évi- 


demment  de  sens  qu’autant 


que  JC  est  plus  grand  que  o,  et  il  vient 


cot.r z= 

.V 


- £ [e  ” + e-«(>  I-')  4-  e-«P 

formule  dans  laquelle  a;  peut  prendre  toutes  les  valeurs  de  - 7:  à + ,1.  Cette 


2. 


dernière  égalité  revient  à 


( 12  ) 


COt  X — - ~ 


f e~^^da\e  e e' ^ J 

Jo 

r e-’d«.  [ 

t/o 


g-K7r_^g-2a;r_,_g-3«:r_^  . . .]; 


OU  bien  à la  suivante 


COt 


x--==f 

./O 


I — e 


-I 


00  ^ — ax  „ax 


dcf.  : 


changeant  j:*  en  - e,  on  aura 


1 2 

COt  - V 

2 P 


ap  «P 

e ^ — e 2 


du. 


Changeons  ensuite  u en  2 x,  du  devient  idx.,  et  les  limites  restant  les  mêmes, 
il  vient 


I 

V 


- COt - e = 
2 2 


r/x, 


formule  qui  subsiste  pour  toutes  les  valeurs  de  e comprises  entre  — 2 t: 
et  + 2 7T. 


7.  Il  est  aisé  de  déduire  de  là  une  nouvelle  expression  des  nombres  Ber- 
noulliens.  En  effet, 


/ p.r 

( — + 

\ I 1.2. 


,3  i.2.3.4.5y 

et,  par  conséquent,  en  vertu  de  la  formule  (6), 

P .1  2 p^  xrfx  2 p‘ 

I — - COt-  V = ~ f 1 / h . . . 

2 2 I Jo  e27TX_j  1.2.3  J^ 

2P’" 

I .2.3.  . .2«  — I Jo  — l' 


( i3  ) 

nous  avons  aussi,  en  vertu  de  la  formule  (,)  où  l’on  change 

i--cot%  = 5lü’  + 5;,  *■ 

O rk  _ _ I / 


e X en  - v, 
2 ’ 


2 2 


'■2  4 1.2.3  6 ■ 771.3.4.5 

Identifiant  ces  deux  développements,  il  vient 


R.= 


xd.T. 


= a.a  r 

Jo 

b,3=2.4  r 
Jo 

H3=2.6  r 
Jo 


Jc^dx 

„'ÎTCX 


2«  1.2 


3.  . .(2  « — j]* 


B„  = 


^ . 2n 


f 

«/o 


=0 


^în— I 


Prenons  maintenant  l’expression 

t/x 

Jo 

dans  laquelle  on  fait  2 tt  .r  = j,  on  aura 


— ~ et  dx 

, 2 TT 


27T 


do. JC  l’intégrale  précédente  devient 


Mais 


par  conséquent, 


j__  J’"-  dj  _ I px. 

T:z7=r-- 


I _ e-7  — i + e^  + e2/4-  e-^r  + . . . ^ 


2^  r /ISO 


'(i4) 

En  nous  reportant  maintenant  aux  propriétés  de  la  fonction  T [intégrale  en- 
lérienne  de  seconde  espèce  (■'“)],  nous  avons 


e~^  Y[nn)—  i.  2.  3 [in  — i ) , 

= — r ( 2 n)  = 1 . 2.  3 (^2«- 

= 1-2.  3 ....(2  7^  - 1 


Donc  enfin 


B„  = 


1.2.3 


]■ 


Telle  est  l’expression  générale  des  nombres  Bernoulliens  développée  en 
série  convergente.  Cette  expression  contient  la  transcendante  tt  et,  en  outre, 
la  somme  d’une  série  indéfinie,  mais  convergente. 

La  valeur  B„  montre  clairement  que  si  n augmente  indéfiniment,  B„  de- 
vient plus  grand  que  tout  nombre  donné. 

Elle  montre  aussi  que  la  somme 


I I I 

^ 


ne  contient,  si  ?i  est  entier,  que  la  seule  transcendante  n,  puisque  nous  avons 
trouvé  une  expression  rationnelle  de  B„. 

Cette  même  forme  fait  connaître  (ce  que  nous  savions  déjà)  entre  quelles 
limites  x doit  varier  pour  que  la  formule  (2  ),  qui  donne  la  valeur  de  tang  x, 

soit  convergente.  La  limite  du  rapport  de  deux  termes  consécutifs  est  — ^5 


et  ce  rapport  sera  < i tant  que  x sera  compris  entre  — ^ 

Ainsi,  bien  que  B„  augmente  indéfiniment,  la  série  (2)  sera  toujour.s  con- 

vereente,  si  x varie  entre  — - et  -t-  -• 

” ’ 22 (*) 


(*)  A la  fin  de  ce  travail  je  traiterai,  dans  un  appendice,  des  diverses  propriétés  des  inté 
grales  eulériennes  dont  j’aurai  eu  à me  servir. 


( i5  ) 

On  verra  de  même  que  la  formule  (i),  qui  donne  la  valeur  de  cot  a.,  est 
convergente  pour  les  valeurs  de  x comprises  entre  — et  + tt. 


III. 

8.  Je  vais  maintenant  établir  entre  les  nombres  de  Bernoulli  des  relations 
dont  je  ferai  usage  dans  la  suite  de  ce  travail. 

Ecrivons  les  deux  égalités  connues 

III 

cot  - e = / 1 

et 


* • . I B,  t» 

cot  - c = — 1 ’ — _j_ , 

2 2 ,.2  I. 2. 3. 4^  12.3.4.5.6 


B„ 


1.2.3... 


2« 


Si  dans  Texpression 


O , 


dx 


nous  prenons  an  - i fois  la  dérivée  par  rapport  à e de  la  quantité  sous  le 
signe  j , il  est  aisé  de  voir  que  nous  aurons  les  résultats  suivants  : 


2.  x"‘-‘ 


cette  dernière  expression,  quand  on  y fait  e = o,  devient 

En  prenant  également  les  dérivées  jusqu’à  l’ordre  an  - i inclusivement, 
dans  le  développement  de^-Icot%,  et  en  faisant  e=o,  on  trouvera?:^. 
Ce  qui  conduit  à la  relation  déjà  obtenue  d’une  autre 


maniéré 


f 


dx 

4« 


Si  l’on  pose  maintenant 

I I 2 cos  V 1 

\ = — cosecot-P  = çe, 


a cos  V 


il  viendra 


( ) 


(7) 


I. 


^ fiVX  fi—fX 


■ cos  V.  dx  = ç)  P. 


On  a d’ailleurs  aussi 


mv  — 

' V 


2 cos  V ^ l 

- cot  - 0 + sin  V. 
2 


Les  développements  en  série  de  cos  v et  de  sin  p sont 


P 

cos  fr»  ==  I h — 


1.2  r.  2.3.4  1.2. 3. 4*5. 6 


sin  P = I — 


.2.3  I .2. 3. 4-5 


Donc 


<p  y 


I t' 


2 


2 1 .2.3  3 1 . 2. 3.4.5  41.2.3.4.56.7 

n — I 


,,2  n—i 


n 


— cot- y, 

I .2.3. . . (2rt  — l)  2 


Et  remplaçant  cot -y  par  sa  valeur,  il  vient 


yv  = B,.-+  + ,.,.3.4.5 


_ /Z  — i\ 

n ■ n j ’ 1.2.3.  . . (2«  — 1) 


Rappelons-nous  maintenant  la  formule 


cos  y = 


.d’où 


e'’^  cos  y = 


k(*+/— 1)  _1_  g»(x— v'— 1) 


( »7  ) 

prenons  les  dérivées  jusqu'à  l'ordre  an  - , inclusivement,  nous  aurons 

—diT-  = â + V'-  i)  + (x  - y/IT7) 

é/’e^^cost'  ir,  , , , , 

~'rfp.  - = 2L('^+  \/+i)  (x  - 


d'”-' e"*  cmv  I f / ^ , 

— 2 • K-^'  + V - -f-  (x  - y'ir7)a"-< 


faisant  ensuite  v = o,  il  viendra 


! d'"-'  e’*  cos  v\ 

\ d^>  )o~ 


_ (x  + yCT r-  + (x  _ 


^(,5/1-1 

et,  en  changeant  a?  en  ~ x, 

ur-' e-"- cosi>\  _ - (x  + - (x  - 

i r/o*^ — * / - — ; 


/ 0 


donc 


frf»"- (e^*cosp  — e->*cosp)-|  r , , 

L rfp’"—  Jo  “ ^ {x  — \/ ~\Y"  ' 

et  comme 


(v/^rr)’"-  = ± ■ , 

V — I 


on  peut  écrire 

jprf»" 


' ( c"  — c— * ) cos  P 


1 (ï  + .*^v/^r”~'-(i-xi/=:Tr"- 

Jo  ■ —jTr 


En  différentiant  — i fois  la  fonction  ipe,  on  a évidemment 
fait  i>  = O,  ’ 


apres  avoir 


rf»"-'  yp\  B 

rfp»”-‘ 


n — 1 


0 « n 


3 


( iB  ) 

et,  par  conséquent,  en  vertu  de  l’équation  (7),  dans  laquelle  on  a pris  les 
dérivées,  par  rapport  à P,  jusqu’à  l’ordre  2«  — i,  et  où  l’on  a fait  v — o, 
il  vient 


(i  -\-x\J — i}’"  ‘ — (i  — X y/— i)’"  ' dx B„  ^ n 

Jo 


n n 


En  développant  la  quantité  sous  le  signe  ^ > et  faisant  les  réductions, 
trouvera  sans  peine 


on 


/-[ 


(2/1  — l)x  (2/Z  — l)(2rt  — 2)(2/J  — 3) 

i TX3  ^ 

(2/Z— i)...(2/Z-5)  (2«~i)(2«-2)„2„_3__2«-I 

' O — / tu  — — cX 

I .2.0.4»5  ï .2 

5i. 

« /Z 


] 


Mais  on  a 


Jr"  X*"-'  B 

-777 dx  — 

0 e ^ — 1 4 


et  si  l’on  fait  successivement  n égal  à i,  2,  3,...,  on  arrive  à la  relation  sui- 
vante : 


(8) 


, xB,  l2n — i)(2n  — 2)  (2/t  — 3)  B, 

( 2 « — I ) ^ — 

' ' l 1.2.3  2 

(2/1—  l)...  (2/1  — 5)]^  _j_  (2« — 1)  (2/Z  — 2)  B„_ 

I.2.S.4<5  3~*~'''  1.2 


B„-,  j 

/»  — I ' 


n — I 

n 


9.  Reprenons  les  formules 


£ 


“ x^”-'dx  __  B„ 

'0  — i ~~  P’ 


et 


£ 


e’*  — e~'- 


O c 


.lit 

dx  = COt  - V ; 


P 2 2 


— 1 


( *9  ) 

par  le  changement  de  en  elles  deviennent 


changeant  ensuite  en  Vy  on  a 


i: 


— e' 


dx  = - — cot  P, 

Cl  ' 


et  multipliant  par  cosv, 


r 


^ ")cOS«'  , COSi»  I 

— r/a:  = h sin  v : — 

e‘  — I t'  sin  y 


Si  nous  prenons  encore  in  — \ fois  les  dérivées  par  rapport  à et  que 
nous  fassions  = o,  nous  aurons 


(co%v  I \n 

h sine 1 

^ c sine/ 

h sj~l  — 

dv'"  ' 

d’ailleurs 


cos  V . I 

h sin  $^  = - H- 

V V 


I . V 
I .2 


et 


1.2. 3. 4“*"  1.2. 3. 4. 5. 6 1.2. 3. 4. 5. 6. 7. 8 


• • • > 


I 

^ — = — cosec  V = 

sine 


2 B,  e 

1 .2 


(2  - l) 


262^’ 

I 2.3.4 


(2^-1) 


263^^  (2>  — 1)  ^ ^ ^ 

1.2. 3. 4. 5. 6 r.2.3...8'‘^ 


par  conséquent, 


( 20  ) 


cos  P 

1-  SIUV 

V SIDP 


I .2. 3.4‘5.6' 


.2.3...2W 


i)i:(2«-  i)]> 


d’où 


[ 


/ COSP  . I \ 

^:n-i  { H Sm  P : 

\ P sin  P / 

'<fp«^ 


1=- 


-î-  - l)±(2W-  l)], 

2n  ‘• 


et  aussi 


£ 


(,  1>«-<  — (i  — :c  v/— . 

“v/=7 


dx  2 ra  — I . I 


+ -Î-.ïB„(2“-'  - l). 


2/j  2n 


Développons 


(i+^\/— 0^"  ^ ^ ^ et  le  premier  membre  de 


f=r. 


l’égalité  précédente  prendra  la  forme 


r*  X J 

[in  — i)  I dx — 

^ Jo  — I 

-- 

1 . 2. 3 -4 *5 


2.(2«  — l)  (2W  — 2)  f 2/Z  — 3) 


X» 

J O — I 


dx 


+ - - Jo  I e--. 


Rappelons-nous  maintenant  que 


r^dx  = 

Jo  — ï 


2’".B„ 

■ ) 

2n 


et  en  faisant  successivement  n égal  à i,  2,  3,...,  il  viendra 


9.//  — 1 f2n-l)(2rt— 2)(2rt— 3)  2‘B, 

TX3  4 


1 . 2 


(2rt  — l)(2« — 2)  2^"  ' Bn_l 
1.2  2 n — 2 


±:^-2B„(2»"  ' — i)4- 

2n  2Tl 


ou  encore 


( 21  ) 


J I 4_  ' 2^  B — — 2)(2/»-3)  2'  B, 

t-3  ’ ■ TITs 

± <i^  0(2/Z-2)  2-'->B„_,  _ ^ _ 2 B, 

1.2  2«_2  — — «;  — 

Multipliant  les  deur  membres  par  ^ . ±,  nous  arrivons  à celte  nou- 
velle  relation  entre  les  nombres  Bernoulliens 

_i ^ L_  A_._ 

.2.3.. .(2/î l)  2'“  1.2.3..  .2/2 2 * 2’"-*  I .2 

_J Bg-t a’"  _ , 

1.2  2’  I.2.3...(2«— 2)  ~ ~ “2”>-  I .2.3...2«' 

Remarque.  — Dans  toutes  ces  formules  le  signe  supérieur  répond  à n pair 
et  le  signe  inférieur  à îi  impair. 


(9) 


1.2.3. 


( 22  ) 


DEUXIÈME  PARTIE. 

FORMULES  D’EULER  ET  DE  STIRLING. 


I. 


Formule  d’Euler. 

10.  La  formule  désignée  ordinairement  sous  le  nom  de  formule  d’Euler, 
et  qui  a une  grande  importance  en  analyse,  est  la  suivante  : 


(lo) 


i - j(Fx  - Fæ„)  + ^(F':r  - F'x„) 


B,  A 


1.2. 3.4^^  ^ i.2.3.4.5.6*^^'''^  ••• 


Cette  formule  est  donnée  dans  le  Traité  de  Calcul  différentiel  et  intégral 
de  Lacroix  (tome  III);  mais,  selon  les  habitudes  des  géomètres  de  cette 
époque,  l’auteur  ne  s’occupe  nullement  de  la  condition  de  convergence,  ni 
de  l’expression  du  terme  complémentaire. 

Elle  a été  reprise  et  perfectionnée  par  Poisson,  dans  un  Mémoire  sur  le 
calcul  numérique  des  intégrales  définies. 

Iæ  célèbre  géomètre  donne  dans  son  travail  une  expression  remarquable 
du  reste  dans  un  cas  particulier.  Cette  même  formule  a été  l’objet  des  re- 
cherches de  Jacobi,  qui  arriva  à une  expression  du  terme  complémentaire, 
pour  le  cas  où,  la  différence  x — x^  étant  divisée  en  parties  égales  à 4, 

^ F(3«)  (X  + 2)  et  (x  H-  2) 

restent  tous  deux  positifs  ou  tous  deux  négatifs  à la  fois,  quand  2 varie 
de  o à A. 

Mais  sa  discussion  n’est  pas  générale,  et  il  restait  toujours  à reconnaître 
s’il  n’était  pas  possible  d’évaluer  le  reste  dans  le  cas  où  la  fonction  et  ses 
dérivées  étaient  assujetties  à la  seule  condition  de  demeurer  continues  entre 
les  limites  dans  lesquelles  la  variable  reste  comprise. 


( ^3  ) 

Le  problème  ainsi  posé  a été  résolu  par  M.  Malmstén,  professeur  à Upsal 
{Journal  de  Crelle,  tome  XXXV). 


11.  Au  lieu  de  présenter  immédiatement  la  formule  d’Euler  sons  la 
forme  (lo),  M.  Malmstén  prend  la  suivante; 


(il)  hf  X— X 


B,  A' 

1 .2.3.4 


etc.. 


d’où  l’on  passe  sans  peine  à la  première. 

L avantage  qu’il  y a à employer  la  forme  (ii)  pour  la  discussion  fonda- 
mentale est  évident. 

En  effet,  si  l’on  considère  l’expression  '^Fx,  de  deux  choses  l’nne  : ou 
bien  elle  représente  une  fonction  continue  de  .r,  et  dans  cette  acception  la 
caractéristique  ^ est  tout  simplement  l’inverse  de  A,  comprenant  implicite- 
ment avec  elle  une  fonction  indéterminée  périodique,  de  telle  sorte  que  V j \r 

est  alors  complètement  indéterminée,  et  il  n’y  a pas  lieu  de  chercher  un 
terme  complémentaire;  ou  bien,  et  c’est  le  cas  le  plus  ordinaire,  on  prend  la 
somme  entre  des  limites  déterminées.  Mais  alors  elle  ne  représente  plus 
une  fonction  continue  de  ;c,  de  telle  sorte  qu’en  la  prenant  dans  cette 
acception,  on  restreint  souvent  la  généralité  des  résultats.  C’est  ce  qui  arrive 

par  exemple,  quand  on  part  de  '^Fx  dans  la  recherche  du  développement 
de  IT  (x  -H  1), 

Il  est  inutile  d'ajouter  que  l’emploi  de  la  formule  (1 1)  ne  présente  aucun 
de  ces  inconvénients. 

Avant  de  la  démontrer,  rappelons-nous  d'abord  un  tliéoréme  connu  lela- 
tif  au  développement  de  Taylor. 

Soit  une  fonction /r  qui  reste  continue,  ainsi  que  ses  2n-hî  premières 
dérivées,  quand  on  fait  varier  de  æ;  à x H-  On  a évidemment 

/(07  + ù)  -yx  =J^  J'  {x  -+-  z)  dz  =£f  {x  + h~  z)  dz. 


{ M ) 

Prenant  l’intégrale  indéfinie  J f {x -h  h-z)  dz,  et  intégrant  par  parties,  on  a 
C -h  J f h—z)dz—  zf'{x  -i-k  — z)  {x  -+-  k — z) 

mVt 

+ + A - 2)  + . . . + , .^.3  /«"'(:»:+  A - Z) 

+ / tÆâX) 

En  intégrant  ensuite  entre  les  limites  o et  il  vient 

t\x  + h)  -fx  = Ayi  = hf'x  + ~J"  X + + • • 


I .2.3. . 


et  si  l’on  change  de  nouveau  h — z en  z, 

Ay-  = A/'^+i^/"^+r^ /'-+•■ 

De  cette  égalité  découlent  évidemment  les  suivantes  ; 

A/x  = hf'x  H-  ^/'-^  + • • • + 

Jo  1.2...(2''  — 0-^ 

A/"  X = V'"  X + . . . + i.2.3.''!.(2«-2)-/‘'”’-^ 

r _ f > (x  + Z ) dz, 

Jo  I.2...(2«— 2/ 

x=  - ^ ~~ï~ 


( ^5  ) 

Si  l’on  multiplie  les  deux  membres  de  la  deuxième  égalité  par  A,  h,  les 
deux  membres  de  la  troisième  par  Aj  les  deux  membres  de  dernière 

il  vient,  en  ajoutant, 


^fx  + A,  k^h^^f’x  + AjA^  Ay"'x  + ..  .4-  Aj„_,  A'*"-' 

hf  X + h'fx  ( j + A,)  + h>fx  ( ^ + A, 


A, 


4 1.2.3 


h^n jr-n) 

,h 


+ 


A, 


I.2.3...2«  1.2.3.. .{2/?  — l)  I.2.3...(2« 2 


+ . . . + An 


jf  fi^n^i)[x-uz)dz 


(h~zy 


+ 


A|  h {h  — z)’" 


4 


1.2.3.  ..2«  1.2.3... (2/2 — l)  1.2.3... (2/2  — 2) 

, Aj„_, /i’"-' (/i  — z) 


) 

] 


Profitons  maintenant  de  l’indétermination  des  quantités  A,,  Aj,  A3,...,  Ajn-n 
pour  poser  les  — i équations 


(.3) 


/ I . 

h A.  = O, 

I A,  . 

, „ O H 1~  An  O, 

1.2.3  1.2  * ’ 

I , A,  , A, 


1 .2.3.4  1 .2.3 


1.2 


4- 


A, 


' 1.2. 3. ..2/2  I.2.3....(2/2— l)  J.2.3...(2/2  — 2) 


4“  ...  4-  Ajn— 1 — 


Ces  équations  sont  les  mêmes  que  celles  qui  se  trouvent  dans  Lacroix, 
mais  nous  allons  en  tirer  des  conséquences  bien  plus  étendues. 

Posons,  pour  abréger, 


rs{x,  h)  = hj' X — kjx  — A,  hkj'x  — ...  — Aa„_,  h^^~' 
Il  en  résulte  ' 


zs{x,k)=  — J 


(A  — z)“ 

1 .2. ..2/2 


A,  h 


(A  — z)’ 


1.2. .,(2/2  — i) 
— z) 


] 


4 


Faisons 


( 26  ) 


. (/i — z)’"  kih{h  — — z)’"~’  k-ih'^  [h  — z)’"~* 

I.2...2«”^  I.2...(2rt l)”^I.2...(2rt  — 2)~^I.2...(2«  — 4) 


Aan-i  ( h — z)' 

I . 2 


et 


, k^h^  {h  — AsA*(A  — z)’"-‘  A(,n_i)  A»"-‘  {h  — z) 

",  1.2... (2/2  — 3)  1.2... (2/2  — 5)  I 


Alors 


zs 


{x,h)=z  — J z)dz[(^{h  — z) 'i({h  — z)\. 


En  développant 


on  pourra  écrire 


(f{h  — z)  + 4i(Æ  — z), 


r rïi 


A, 


2/2  1.2. 3. ..(2/2 — l)  1.2. 3. ..(2/2  — 2) 


..(2/2  — 3) 


-h 


kzn — I 


'] 


A^"“'zj  1.2. 3. ..I 

L+ — r 

L 1.2. ..(2/ 


A, 


2/2 — i)  1.2. 3. ..(2/2 — 2)  1.2. ..(2/2 — 3) 


.2. ..(2/2 — 4) 


ç — z) + — ^)  ■ 


r î 

1.2. 3. ..(2/2  — S 
l_  1.2. ..(2/2 — 4j 


-+-...  H ^ h Aj„_, 

A. 


] 


2/2  — 2)  1.2. ..(2/2  — 3) 


A}„_3  . 

• • H :: F Aj„_2 


/i»z2"-3  r I A,  A,  . T 

— ^ 5-:  I 5 H h 1~  A3  I 

1.2. ..(2/2  — 3)1_1.2.3  1.2  1 J 


h‘‘  z’"“ 

1 .2...(  2 n 


T-- -^^'1 

1.2...(2/2—l)  Ll  J 


+ 


I .2.3...  2/2 


!» 


( 27  ) 

etayantégard  aux  égalités  qui  ont  lieu  entre  A,,  ,1  vient,  puisque 

A.  = -i, 

z’"  A, 


I.2.3...2/2  I.2...(a«  — () 

<p(A-z)  + 4/(^-z)=/  + — +...  , 


i.2...[2n  — 2] 
A,  /<’  z’"~’ 


1 .2 

A 5 /i‘  z^"— ‘ 


1.2.3...(2«  — 3)  I.2...(2«  — 5) 


...+ 


A„,_, 


^ /’ 


mais  la  première  partie  du  deuxième  membre  de  cette  égalité  est  évidemment 
9 Z,  la  deuxième  — i^z,  donc 


(p  {h  —■  z)  -h  ^ {k  ~ z)  = (p Z — (il Z, 


et  faisant  z = ~ h. 


d’où 


4'(ï^)=°'  , 

(2^)  nulle  pour  toutes  les  valeurs  entières  de  n,  en  donnant 

à n successivement  les  valeurs  a, '3,  4,  etc.,  on  aura 


et 


donc 


et,  par  suite, 


A3=o,  As  = o,...,  A2^_,  =0, 
<ipz  = <]p[h  — z)  = o; 

9 Z = <p(h  - z). 


ST 


>,//)  = _ f + z)cpzdz. 

«/O 


D ailleurs  pour  trouver  la  relation 


qui  existe  entre  A,,  A3,  A„...,  et  les 

4. 


( ^8  ) 

nombres  Bernoulliens,  il  suffit  de  se  reporter  à la  dernière  des  équations  (r  3) 
qui,  en  vertu  des  relations 


A,  1 A3  A5  . . . Ag„_4  O, 


devient 


A(2n— 2) 


n — I 


1.2  I .2. 3 .4 


1.2, 3.. .(2« 2)  1.2. 3.  ..2/2 


multipliant  par  ar(an),  il  viendra 


n — I 


= aAj(2«  — 1)  + 2k^[in  — 2)(a«  — 3) 

2A2„_a(2  7^  — i)(2n  — 2)...4.3; 


ou  encore 


n — I 


A,  _ (2/2  — 1)(2//  — 2)(2/2 — 3) 
1.2.3 


(an—  i)-  ï -2-  Y . ..  2 I -2.3.4 


A, 

2 


(2/2-i_)^_^^  ,.2.3. ..(an-  2)^. 

^ 1.2. 3. 4.5  ^3  1.2  / n—i 


.3.4.5 

Si  l’on  compare  cette  formule  avec  la  formule  (8),  on  aura  évidemment 

Bs 


* B,  A _ B 

Ao  — — ’ A-4  — — 


1.2 


1.2. 3. 4 


— 


I.2.3.4.5.6"’’ 


(i3') 


A “f* 




B<t— 1 


1.2. 3. ..(2  //  — 2) 

selon  que  n est  pair  ou  impair,  ce  qui  est  la  relation  qu’il  s’agissait  d’établir. 


12.  Considérons  maintenant  l’intégrale  J (pz.dz.  Je  dis  qu’on  peut  écrire 

h 

J'  fzdz  — a l fzdz» 

0 do 


En  effet 


(l>z=  f{h  — z), 


on  a d’ailleurs 


( 29  ) 


£^zciz  = £ fZdz  + £\(h-z)dz. 

a 

Si  l’on  pose  h~  z = z',  \[  vient 

dz=z  — dz', 

et  les  limites  - et  h se  changent  en  ^ et  o;  donc 

h 

jT  (f{h  — z)(iz=.  J'  ^z'dz', 

2 

et,  par  conséquent, 

h 

jf  f2rfz  = 

^zdz. 

La  valeur  de  <pz  est  donnée  par  la  formule 


<pz  = 


1.2.3...2/Î 


A, 


1.2.3.. 772.3...(2«-2) 


A,,_,  ^«-*2= 


1.2 


et  s.  nous  intégrons  entre  les  limites  o et  h,  il  viendra 

f <pzdz  = h^’‘->-il  ^ ^ A,  A,  AT 

' L 1-2.3.. .{2/1  + x)  I.2.3.. . 2/1  77X7(77:^7) 

1 intégration  entre  les  limites  o et  donne 

i f'  (pzdz=/d"^^\  1 -..JL,  A,  I A,„_,  il 

Jo  Li.2.3...(2/H-i)  2='>~*~  I.2.3...2/i'7^“‘'---~*"  T • 

î ai  J 


Dans  la  dernière  des  formules  (.3)  changeons  « en  « + i et  rappelons- 
que 


nous 


A3=  A, 


Aj;,—  , O, 


on  aura 


( 3o  ) 


I.2.3...(2/2  -l-l)  1,2. 3... 2/2 


ri:3 


Aon  O . 


Mais  on  a aussi 


j(  cpzdz  = 2 I (pzdz, 
O J O 


d’où  résulte  la  relation 


(i4)  — Aa„  = 


1 ! _ 

.2.3...(2/2  + i)  2’"  1.2.3...2«  2 


I 

1.2.3  2' 


Propriétés  fondamentales  de  la  fonction  <pz. 

15.  I®.  La  fonction  ç)  Z ne  change  pas  de  signe  entre  les  limites  z=o  et  z=^h 
de  l’intégration,  elle  est  toujours  positive  entre  ces  limites  si  n est  pair,  et 
négative  si  n est  impair. 

Pour  le  démontrer,  remarquons  d’abord,  comme  cela  a été  dit  plus  haut, 
qu’on  obtiendra  toutes  les  valeurs  que  peut  prendre  la  fonction  <pz  entre  o 

et  h en  faisant  varier  seulement  z de  o à - â. 

2 

Si  nous  considérons  maintenant  l’expression  z -H  A,  Æ , dans  laquelle 
A,  = — il  est  évident  qu’elle  sera  négative  pour  toutes  les  valeurs  comprises 

entre  o et  - Æ. 

2 

Mais  pour  « = i on  a 

z’  Al  hz 

(fZ  — 1 

'1.21 

Sa  dérivée  est  z -{-  A,  et  négative  pour  toutes  les  valeurs  de  z com- 
prises entre  o et  ^ ^ ; donc  <p  z décroît,  et  comme  elle  est  nulle  pour  z = o, 
elle  sera  négative  pour  les  valeurs  de  z considérées. 


Si  n = 


( 3i  ) 


2, 


(pz  = 


A,h^z' 


1.2.34  1.2.3 


+ 


1.2 


et 


A,hz^  . 


Mais  l’expression  î + ^ est  négative  pour  toutes  les  valeurs  de 
entre  o et  et  il  en  sera  de  même  de  p'  (i  + ^era  la . 

d’un  nombre  infini  de  termes  négatifs  infiniment  petits.  Or 


Z comprises 


somme 


on  a 


ce  qui,  en  vertu  de  la  formule  (i/f)?  devient 


f — 

Li.2. 


A,  hz 

3 1.2 


Aj  h 


et  comme  l’intégrale  définie  est  essentiellement  négative,  la  quantité  entre 
parenthèses  est  positive. 

Ainsi  qu  il  est  aisé  de  le  voir,  cette  quantité  n’est  autre  chose  que  -f'z; 

la  fonction  ç'2  est  donc  positive  pour  toutes  les  valeurs  de  z comprises  entre 

O et  Comme  d’ailleurs  çzest  nulle  pour  2 = 0,  la  valeur  qui  répond  à 

Al  = 2 sera  positive  pour  toutes  les  valeurs  de  2 considérées,  puisque  la  fonc- 
tion est  croissante. 

Si  Al  = 3, 


et 


( 3a  ) 


9'z  = 


1.2.3.4'S 


Ai^z* 

1 .2.3.4 


A, 

1.2.3 


“4“  h Z. 


Reprenons  l'expression  ^ + AjA’z,  qu.  est  positive  pour  toutes 

les  valeurs  de  z comprises  entre  o et  ^ A.  Multiplions  par  et  il  est  évident 


f-i'lix 


3 1.2 


sera  positive. 

Cette  intégrale  définie  devient 


Z»  Al  Z hr'  ^ ^ . 

ïT^Ts  1 .2  4 * ^ ’ 


donc 


1 A|  z’  A 

1.2.3  5 1.2  4 


Ajz  ^ 

I 3 


est  positif  pour  toutes  les  valeurs  de  z comprises  entre  o et  i /*. 
Il  en  sera  de  même  de 


f z’  1 A,  z'  h A,  z , 


qui  est  évidemment 

/ z‘  A,Az’  AjA^z’X  7«  / 1 * I 

“ 1.2.3  Vi.2.3.4.5  2‘  1.2.3.4  2 

On  pourra  donc  écrire,  en  vertu  delà  relation  (i4), 


2’  I .2.3  2’ 


À\  z>  1 AlZ’A  A,zA’\  . 


( 33  ) 

et  comme  la  valeur  de  l’intégrale  est  essentiellement  positive,  la  quantité  entre 

parenthèses  doit  être  négative.  Mais  cette  quantité  n’est  autre  chose que-y'z; 

donc  ô Z est  négative:  or  pour  z = o,^z  = o;  donc  çz,  fonction  décrois- 
sante, est  négative  pour  les  valeurs  de  z considérées  lorsque  w = 3. 

On  poursuivra  sans  peine  le  même  raisonnement  jjour  toutes  les  valeurs 
entières  de  /i. 

La  première  propriété  énoncée  se  trouve  ainsi  démontrée. 

2°.  La  fonction  fz  a entre  z = o et  z = A un  seul  maximum  pour  z=-/isi 

n est  pair,  et  un  seul  minimum  pour  les  mêmes  valeurs  de  z si  n est  impair. 
Cette  propriété  se  démontre  très-facilement,  car  on  a 


cpz=cp(h  ~ z), 


et  en  prenant  les  dérivées  cp'  z = — (p'(^k  — z], 


d’ailleurs  pour  z = -h^ 


Mais,  toutes  les  fois  que  n est  pair,  ip'z  est  positive  pour  les  valeurs  de  z 
comprises  entre  O et  ^ ^ qui  est  positive  et  croissante  pour  de  pa- 


reilles valeurs  de  z,  croîtra  donc  de  z = o à z = 1 /^,  et  comme  la  fonction 

repasse  par  les  mêmes  valeurs  pour  des  valeurs  de  z équidistantes  de-// 

. 2 ’ 


elle  aura  un  seul  maximum  qui  arrivera  pourz~-/i. 

Un  raisonnement  identique  avec  le  précédent  montrera  que  la  fonction  ^z, 
toujours  négative  dans  le  cas  de  n impair,  si  z varie  entre  o et  h,  passera  par  un 

minimum  pour  z = -^,  minimum  qui  répondra  d’ailleurs  à un  maximum  en 
valeur  absolue. 


14.  Ces  propriétés  de  la  fonction  ^z  une  fois  établies,  revenons  à la  for- 
mule (12). 


0 


Rappelons-nous  que  l’on  a 


( 34  ) 


A3  — A5  — — ...  Ao«_(  — O, 


A,  — 5 A 2 — — = 

2 ■‘1.2 


B.-  . ^ B3 

1.2. 3. 4 ® 1.2. 3. 4. 5. 6 


etc. 


Rappelons-nous  aussi  que  l’intégrale  définie  qui  entre  dans  le  second 
membre  peut,  d’après  ce  qui  a été  dit  plus  haut,  être  mise  sous  la  forme 


/•ft 

j J + z).(p z.(h, 


et  alors  la  formule  en  question  devient 


(■5) 


1.2.3...  1 2 


n— 2 , /»  A , 


Mais  <pz  garde  le  même  signe  pour  toutes  les  valeurs  de  z comprises  entre  o 
et  h',  donc  l’intégrale  définie  est  égale  à la  somme  des  facteurs  (fz.dz,  mulli- 
pliée  par  une  valeur  moyenne  comprise  entre  la  plus  grande  et  la  plus  petite 
valeur  que  peut  prendre  la  fonction  [x  z)  quand  z varie  de  o à /?, 

cette  fonction  étant  continue  entre  les  limites  de  l’intégrale.  On  pourra  d’après 
cela  écrire 


r J + z)  dz  = y ( JT -+- ô/?)  r ozdz^ 

Jn  Jo 


Q étant  compris  entre  o et  1 , et  en  désignant  par  R le  reste  de  la  série, 
R = — y (x -t- ô//)  J (ozdz. 


D’ailleurs 


zizdz  = 


I 

\I  .2.3... 

f 2 « -t-  I j 

“t” 


1.2.3...  2 « 


A, 

: .2.3...  ( 2«  — I ' 


A,  n_A 

r^.3/ 


J 


( 35  ) 


d’où,  ei)  ayant  recours  aux  équations  (i3)  et  (i3'), 


/•  h 

I <p  zdz  = — A2,,  = rt . 

•/ O 1.2.3...?!// 


,2/!+I  . 


do 


ne 


+ SA), 

et  par  conséquent  aussi 


(.6) 


j hJ'.x=àJx-'-h\J\r+ 

I 


R . 

1.2.3...  (2  /i  — 2)  ^ ^ 

B 


Si  dans  la  formule  (i5)  nous  faisons  x = x a\.  jc  = Xo  et  que  nous  retran- 
chions membre  à membre  les  équations  obtenues,  il  viendra 

i A(/'^'-y'.r,)=A>-A>,-A/,(A/a:-A/'a-„)+^’(Ay'.r-A/X) 
('7)1  - - .)  ( 


et  comme  çpz  garde  le  même  signe  pour  toutes  les  valeurs  de  z comprises 
entre  O et  h,  le  dernier  terme  du  second  membre  sera  égal  à fzdz  muiti- 

t/o 

plié  par  une  valeur  moyenne  que  prendra  l’expression 


~f^^n^'\{Xo+  z), 

lorsqu’on  y fera  varier  z de  z = o à.  z = h;  on  aura  donc 
/ \fjc—âijx^-  ^h{Af'.x~-A f'x^] 


B,// 


^{Arx-Af'x„] 


B„  4'"+' 


772:3TT^  + Oh)  - + 5/0]. 


('7') 


( 36  ) 

Pour  arriver  à la  formule  connue  sous  le  nom  de  formule  d Euler,  il 
suffit  de  poser 

+ — z), 


d’ 


ou 


J'  (x  -h  z)  = (x+  z)  -+-  C, 


et  par  suite 

/'  (x  + z)  — /'  (Xo  + z)  = 


c disparaissant  si  l’on  suppose  x — x„  égal  à un  multiple  exact  de  h : on 
a aussi 


fx 

Fxc/x, 

A/'x  — A f’Xo  = F X - F Xo, 
Af"x  - A f"  x„  = F'x  — F'Xo, 


et 

^ + z)  — + z)  = 2]  "*■ 

Donc 


i /?  V Fx=  r F xdx (Fx  — FXq)  + — - (F  X Fx^)... 

i t r>'\  J ^ 

1 T>  /,2n— 2 ,,  X B„  /i’"”*"'  T’2H/  ùî\ 

+ .(F(2«-3)  ^ _ F<2”-3)  xjzp  , " -y  f h) . 

1.7.. ..{ïn — 2)^  JC„ 


Ce  qui  est  la  formule  d’Euler  avec  son  terme  complémentaire. 

Posons  X — x^=ih.,  et  désignons  par  M la  plus  grande  valeur  numérique 
de  F<-">x  quand  x varie  de  x,,  à Xp  + ih,  il  est  évident  que 


y""  F<“")(x  -F  Qh) 


( ^7  ) 

sera  en  valeur  absolue  moindre  que  M^et  par  conséquent 
il  vient  alors 


— tt: ~[x  — xJ  M, 


1.2...  2n 


R)  désignant  la  valeur  absolue  du  reste. 


lo  Supposons  ma,ntcn.int  que + j)  conserve  le  même  signe 
quand  on  fait  varier  z de  o à /j,  ^ 

Dans  cette  hypolhèse/>-„(^+,,  y, 

facteurs  tous  de  même  signe  + multipliée  par'  une  valeur 

moyenne  que  prendra  la  fonction  yz  quand  z variera  de  o à /<,  on  à 

? (ÔA)  X + z)  ifz  = A)  (n:+A)  — yn«) 


donc  en  revenant  à la  formule  { , 5)  on  aura  pour  l’expression  du  reste 

R z=  — ç)  (0^)  ^ y-2«  . 

mais  la  plus  grande  valeur  de  la  fonction  cp  [z]  est  ? donc 

R=—  y (2«) 

9 étant  compris  entre  o et  i. 


La  formule  (i3)  qui  donne  <pz  devient  pour  z = - h 

2 ) Li-2.3...2«  2^’'^  1.2.3...(2n-l)  2^n-, 


( 38  ) 


ou  avec  les  nombres  Bernoulliens 


I I 

1.2.3...  2 fl  2"* 


2 1.2... (2« — i)  2’"~'  1.2  1.2... (an  — 2)2'"“^ 

_ B, I I 

1.2. 3. 4 1.2.3... (2«  — 4) 

B„-, L.-b- 

1.2.3...  (2  « — 2)  1.2  2’’ 


donc  en  vertu  de  la  formule  (9) 


2'''‘  — 1 
1.2...  2 n 2’"’^' 


et  en  se  reportant  à la  valeur  de  R, 


R = qr  0 


1.2...  2 rt 


A j X , 


ce  qui  est  une  nouvelle  expression  du  terme  complémentaire  de  la  for- 
mule (i5). 

De  même  si  dans  la  formule  (17)  JV*»  {x  + z)  -/  = ~'(x„4-z)  con- 
serve le  même  signe  quand  on  fait  varier  z entre  les  limites  o et  h,  on  a 

z)  — -I-  z)]<pzr/z 

= (-^0  + /i)  X +y-^"xo] 


D’où  résulte  l’égalité 

[ h{fx -j'x,) = Ayx  - Ayxo-  ^[a/'x  - A/'xoj  -t-  ^ [ù.y  "X  - Ay".v„J 


(18)^ 


r ^ r-(2„_2)  ^ / i2«-  2)  J 

1.2.3...W— 2L  J 

Uy(2«)^_  A/(=*")xoj. 


G 2’''  — 1 


B„  h^-' 

1.2  . . . 2 n L 


10.  Pour  repasser  à la  formule  d’Euler,  on  fera  comme  au  para 


graphe  14-,  et  il  viendra 


( % ) 


>9)' 


'X  ; (f-  - F-«)  + ^ (r  .r  - P X.)  - 

R /j-'n— î / 

± r 

I . 2 . . . . ( 2«  — 2 y '^o 

^ B„  h'-''  2"'  — I / 

-I-  ô ( p'(2"-n  p(2/f-)) 

I 2...7..’l  2-"-'  \ -L  I .tj 


Il  est  presque  mutile  d’ajouter  qu’eu  employant  celte  formule  on  suppose 

que  h (2„)  conserve  le  même  signe  quand  on  fait  varier  z de  z = o 

■à  Z ~ h. 


Je  fera.  ,en,a,q.,er  .ci  q..e  po...-  a.river  a..  ler.,.e  co.npléme.,tai.-e  de  la 
sene(io')  a..ss.  b.ei.  q..a  ceh.i  de  la  série  (.9),  il  n’est  pas  nécessaire  de 
supposer  , ent.er.  Seule.., ent,  dans  le  cas  of.  i se.-ait  fractio..nai.e,  chac... 
de  ces  restes  devrait  être  accompagné  d’nne  quantité  C=Ç.r— Ç.r  étant 

une  fonction  péiiodique  de  a:  q..i  repasse  par  les  mêmes  valeurs’q..a..d  0.. 
change  jc  en  x h. 


Nous  n'avons  parlé  de  ce  cas  que  dans  le  but  de  mo.,t.-er  la  généralité  ,1e 
la  méthode,  nous  le  laisserons  de  côté  désormais. 


17  Poisson  trouve  un  terme  complémentai.e  dont  la  forme  diffère  nn  peu 
de  cel  e qn,  précédé,  et  il  arrive  à cette  conséquence  remarquable,  savoir  que 
s.  F X garde  le  meme  signe  depuis  x = x.  jusqu’à  x = x,  le  .-este  de  la 
serie,  en  s arrêtant  a un  certain  ternie,  est  moindre  que  le  dernier  conservé. 

On  ai  rive  ici  a la  meme  conclusion. 

En  effet,  revenons  à la  formule  ( , 5),  dont  le  dernier  terme  est 


%/0 


Si  nous  supposons  que/«-)(x  + Z)  garde  le  même  signe  entre  les  limites  de 
I intégration,  il  sera 


R = qzô. 


1.2...  2.n 


X,. 


( 4o  ) 

(-t  l’on  pourra  mettre  (i5)  sous  la  forme 


/•  h . El  ^1*’  A /V/  -t-  B„_|  . /(2n-2 

h.f'x  — ùijx  — -A/  o:+  — A./  ac:  j_2....(2«  — 2) 


(20) 


B„  /«• 


1.2...  2 « 


■(2«) 


. X. 


d’où  l’on  tirera,  comme  pins  haut, 

h if'x  -f'Xo)  = A/x  - ^fx^  + . . . , etc. 

et  enûn  ^ 

hY^x=^  r Yx.dx-l[Yx-^x^)-  ■ • 

/p(2/i-i)  y.  _ p(2«-n 

1 . 2 . . . 2 « ' 

^ /'i_g  -Ai^(F(2'^-'>x-F(=‘"-').r)o, 

— y ^ 2^"-'  J 1 .2 . . . 2 rt  ^ 

où  l’on  suppose  que  (ar  + 2)  ne  change  pas  de  signe  quand  ^ varie 

de  Z = O àz  = ù. 

Cela  arrivera  évidemment  si  F x garde  toujours  le  même  signe  quand  x 
varie  de  a,’  = a:o  à x = x,  comme  le  suppose  Poisson,  mais  pourra  encore 

arriver  dans  d’autres  cas. 

Si  nous  remarquons  maintenant  que  est  plus  petit  que  2,  1 — & 2-"-'  ' 

sera  en  valeur  absolue  plus  petit  que  i,  et  le  reste  de  la  série  sera  évidem- 

ment moindre  que  le  dernier  terme  conservé,  sans  qu’il  soit  possible  par 
cette  analyse  d’en  déterminer  le  signe. 


{ 4i  ) 

^ 18.  Cela  posé,  M.  Malmstén  se  place  dans  l’hypothèse  (æ'+  z)  et 

f [x-hz)  ne  changent  pas  de  signe  quand  z varie  entre  o et  h.  Il  peut 
alors  se  présenter  deux  cas  : ces  deux  fonctions  ont  le  même  signe,  ou  elles 
ont  des  signes  contraires;  et  dans  chacun  de  ces  cas  il  trouve  une  expres- 
sion du  terme  complémentaire  avec  son  signe  en  évidence.  Reprenons  d’abord 
la  formule  (i6),  et  changeons-y  n en  n-hi. 

Le  terme  complémentaire  deviendra 


1.2.3...(2«-}-2) 


/ 


2/Î4-3 


d/l). 


Ce  terme  dans  la  formule  (20)  est 


, 2="  — I \ „ 


car  on  sait  que 


■(2/J+U 


(x-hz)ffz= 


mais  il  est  évident  que  si  (a:  + z)  et  (x  + z)  gardent  le 

même  signe  quand  2 varie  de  o à i - 0 sera  positif,  soit 


ô. 


2-«—  I 
2’"-' 


0,; 


et  dès  lors 


_ 9,  B„  //’" 

I . 2 ...  2« 


A/OT(x) 


représentera  le  terme  complémentaire  de  la  série  (i5). 

§i  jT{2n+3)  2]  gjy(2«+o  sont  de  signe  contraire  entre  les  mêmes 

limites  de  2,  I — 0.  ^ sera  négatif  et  égal  à 


-02. 


2*"  — I 
2^''-' 


6 


(40 

S-i  étant  compris  entre  o et  i,  et  le  reste  de  la  série  mise  sous  la  forme  {20) 


sera 


i .2.3  ...  2 n 2 
on  pourra  donc,  selon  les  cas,  écrire 

( , 5' ) 4/-X  = -1/—  ; A/' X + . . . . ± B„,  ) A/-' + 6,  A/» (x) 


OU 


(20')  hf'.v  = Afx ùif'x 

2 


R Ain  0.,  R /)2K  o2«  I 

■A/^"x=p--  2 iA/».r. 

■ ^ 1.2.3  ...  2/2  2’"~' 


1.2...  2 n 


Revenons  maintenant  à la  formule  d’Euler,  et  dans  (10')  changeons 
n en  n+  i,  le  terme  complémentaire  deviendra 


/2’"+’  B„ 


1.2.3  ...(2/24-2) 
et  dans  {21)  ce  même  terme  est  de  la  forme 
B„  A’» 


^ + Q h). 


1 . 2 ...  2/2 


-«■O)  X' 2?'“’ 


car 


Ç ^ p(2«)  (x  4-  z)  X - 

Jo 

En  raisonnant  comme  plus  haut  on  verra  que  si  ^ F*”'*'*  (x  4- z)  et 

2 F [x  4-  z)  gardent  le  même  signe  pour  toutes  les  valeurs  de  z comprises 
entre  o et  A , on  pourra  écrire  les  deux  formules  suivantes,  où  les  signes  des 
termes  complémentaires  sont  connus  : 

1°.  Si  (x  4- z)  et  51  F*"(x-t-z)  sont  de  même  signe, 

Fo:=  Ç Fx.<ir  — - (F-c  — Fxo)4-— ^ — (F'x — F'x<,) ^(F"'x  — F"'j:o) 

X 2^  ‘1.2^  1.2.34' 


(22) 


B„_,  /2^"-» 

' 1 .2  ...  (2  n — 2) 


6 R /i*'* 

(F  X—  F Xo)  ZC ^ ( FO"-'  ) X — FO"-'  ) X J. 

^ ^1.2. ..2/2 


2' 
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».  s.  et  2' F'“'(^+z)sont  de  signes  contraires, 


I .2.3.4 


y \ — -, fF^"-3  r — F’"-3 

1 I.2...(2«  — 2)  ^ 

f _ B„,  /i^- 


1.2  ...  2n 


(F^-'j:— F("‘-‘).r„):5IÔ  — i B„4> 


I .2...  2« 


Dans  (aa)  et  (23),  on  suppose  la  continuité  de  la  fonction  F(ic+  a)  et 
e ses  denvees,  continuité  qui  doit  s’étendre  jusqu’à  F {x  + z] 

(-es  formules  résument  le  perfectionnement  remarquable,  ou,  pour  mieux 
d'EulIr  ‘‘«fi"*''''  <I“>^  M.  Malrastén  a donné  à la  formule 

Les  séries  trouvées  plus  haut  peuvent,  dans  un  grand  nombre  de  cas  pré- 
semer  au  commencement  une  convergence  sensible,  puis  devenir  ensuite 
ivergentes  par  suite  de  l’accroissement  des  nombres  Bernoulliens,  mais 
expression  SI  simple  qu’on  a donnée  du  terme  complémentaire  permettra 
en  general  d assigner  dans  tin  cas  donné  le  terme  auquel  on  devra  s’arrêter 
pour  avoir  la  plus  grande  approximation  possible. 


IL 

Formule  de  Stirling. 

D 

19.  Parmi  les  nombreux  développements  en  série  que  donne  Stirling  dans 
son  ouvrage  intitulé  : Methodu,  dijjerentialis , sine  Tmclatus  de  sumniatione 


( ) La  formule  (22)  avait  été  trouvée  par  Jacobi, 
M.  Malmstén. 


mais  la  formule  (23)  est  due  à 


6. 
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serierum,  se  trouve  la  suivante  : 


(^4) 


1.2  X 


B 


lin  + xlx  — X -\~  -Ix 

I ^ I 

I X 3.4  5.6  x^ 


A.  ‘ 

5.6  x^ 


Cette  formule,  qui  a gardé  le  nom  de  son  auteur,  est  justement  célèbre,  tant 
par  son  utilité  que  par  la  propriété  singulière  qu  elle  présente. 

Elle  donne,  en  effet,  le  moyen  d’évaluer  rapidement  la  somme  des  loga- 
rithmes de  la  suite  des  nombres  naturels  depuis  i jusqu’à  x,  et  avec  d’au- 
tant plus  d’exactitude  que  ^ est  plus  grand.  Cette  opération  se  présente  sou- 
vent dans  le  calcul  des  probabilités,  et  il  serait  la  plupart  du  temps  impossible 
de  la  faire  directement. 

Mais  ce  qui  a surtout  attiré  l’attention  des  géomètres,  c’est  que  cette  série, 
très-convergente  pour  de  grandes  valeurs  de  x quand  on  se  borne  aux  pre- 
miers termes,  devient  divergente  quand  on  en  prend  un  plus  grand  nombre, 
et  cela  quelque  grandes  que  soient  les  valeurs  attribuées  à la  variable  x. 

Malgré  cette  singularité,  elle  n’a  pas  été  abandonnée  des  analystes;  ils  1 ont 
conservée  à cause  de  sa  grande  utilité,  et  en  ont  légitimé  l’emploi  en  donnant 
son  terme  complémentaire. 

C’est  ce  qu’ont  fait  entre  autres  MM.  Liouville  et  Cauchy. 

Nous  allons  traiter  la  question  en  donnant  d’une  manière  générale  le 
développement  de  i),  où  x peut  avoir  des  valeurs  positives  quel- 


Ce  développement  est  une  application  de  la  formule  (i5')  trouvée  plus 


conques. 


haut. 


20.  Pour  y arriver,  partons  des  formules  connues 


et 


Jî 
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fl  après  cela,  nous  pouvons  écrire 


J'x=l.T{x)=z  r“  [(or  - i)e-='  -îlllllîZlfl 

4/0  L I — e~~“  * I 

et,  en  prenant  les  dérivées  successives  par  rapport  à x. 


_ d.ir(. 


J"x=- 


<lx 


J'"' X r^a.f—ax  ^ 

J O t — e—« 


OC  ) 

r 

0 i—e- 


cl  a. 


~ — r” 

Jo 


• f/a. 


selon 


que  i est  pair  ou  impair. 
Cela  posé,  on  a 


Ay^r — J ' {x z).dz  =.  J'  l.T{x+z)dz, 

et  si  l'on  fait  h = \ , 


A Jx  = f‘  /.r{x-t-z). 

' Jo 


dz. 


Posons 


d’où 


X+Z^J 


dz  = dy. 

Les  limites  de  l’intégrale  deviennent 


d<x 

> 

(X 


X — I 


et  x; 
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et,  par  conséquent, 

I.T{J -^\)dy  = l.T[f + \)dj  - l.T{r-^\)dy 

= C l .r{j  + i)dy  + C l.T{jr-h\)d/  - r l.T{jr  + i)dj. 
J O 

Soit 

y ^ I = y,  d’où  dj  = dj'; 
les  limites  o et  x — i se  changent  en  i et  x,  et 

jf  + = i.T{y).dy-, 

supprimant  l’accent,  il  vient 

A/or=  jT'  i.r(jr+l)rf/+  J [l.T{jr+,)-  lTf]df 


mais 


r.(j  + 

r(j) 


= J, 


et,  par  conséquent, 

àfx=  f l.T  {j i)  dj -h  r ly.dj=  Ç l.T[jr-^i).djr->rxlx—x+\ 

J O •-'I 

D’après  une  formule  connue,  on  a 
Jo 


(a) 

Donc 


i^fx  = - Z.  2 71  -+-  xlx  — a 


et  en  prenant  les  dérivées  successives, 


A/'".r 

A 


~ Ix^ 


I . 2 


A f^’^x 


I .2.3.4 

: î 


^fl.in-2)  _ 

A/=*”.r  = 


1 .2.3..  . (2  « — 4) 

3 

1 .2.3..  . (2 /?  — 2) 


Revenons  à la  formule  (i5'),  et  rappelons-uous  que 2)  et 
y(2«-v3)^^  2)  sont  essentiellement  positifs.  Il  vient,  en  faisant  h=n^ 

/ 'O-  = Ayi-  _ i A/-^  + A 


J .2.  . . (2/î  — 2) 


1.2. 


2/2 


A./‘»x, 


et  remplaçant,  dans  la  relation  précédente,  les  A par  leurs  valenrs,  nous 
trouvons 


= ; /27T  + (x  - i)  /x  - X -h  A I 

\ 2/  I .2  X 

Remarquons  maintenant  que 

+ i)  = a:r(.r), 

d’où  ' 

l.T{x  4- i)  = /r(x)  + /a-. 


et  nous  aurons 


lT{x 


\)  = '-l2n  Ix  -^-\lx  — X -h  ~ 

^ \ 2/  1.2  X 3.4  -r’ 


_ e.B„ 


(2/2  3)(2rt — 2)  x’"  3 ”*■  (2/2  — 1).2/2  .r’"— 


Il.1 

5.6  x‘ 


Ainsi  que  nous  l’avons  dit  plus  haut,  la  formule  de  Stirling  est  rapidement 
convergente  dans  ses  premiers  termes  pour  de  grandes  valeurs  de  x,  mais,  en 
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s’avançant  plus  loin  dans  la  série,  les  termes  ne  tardent  pas  à croître  de 

nouveau,  et  cela  quelque  grand  que  soit  x. 

Pour  s’en  convaincre,  il  suffit  de  prendre  le  rapport  d’un  terme  au  prece- 
dent et  de  passer  à la  limite  en  faisant  n infini. 

Ce  rapport  est 

3)(2«  — 2)  1 

“ B"~'  (2/2  — l)  2rt 

qui  doit  à la  limite  être  plus  petit  que  i,  ce  qui  conduit  à l’inégalité 

_ (2«— 2)(2/J-3) 




ou  à une  valeur  infinie  de  x quand  n est  infini. 

Néanmoins,  puisque  l’on  connaît  le  terme  complémentaire  de  la  série,  il  sera 
toujours  facile,  dans  un  cas  donné,  de  trouver  le  terme  auquel  on  devra  s ar- 
rêter pour  avoir  la  plus  grande  approximation  possible. 

1 


Développement  de  Lt[x)  en  série  convergente. 
21.  Nous  partirons  encore  de  la  formule  connue 

P Q — a - — e — « ^ n d<x. 

qui  donne  en  prenant  la  dérivée  par  rapport  à x, 

r'(^)  ^ f"  flZl  _ r/a. 

Jo  L « i-e  «J 


Posons 

d’où 


e-“  = t, 

da  — — jdt  et  a—  l.t, 
pour  a = O,  < = ) et  pour  a = co  , i = o. 
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Les  limites  de  l’intégrale  deviennent  donc 


f et  O,  et,  par  conséquent, 


dt. 


Si  nous  faisons  maintenant  = t,  il  viendra 


= -C, 


C étant  une  constante  dont  la  valeur  a été  calculée  par  Legendre  et  dont  il 
suffit  ici  de  constater  l’existence.  Nous  aurons  d’apres  cela 


Or 


iTf) 

r(^) 


-t-  C = 


I — 1^~' 


1 — t 


dt. 


I 


I—  r 


= I -h  -h  t^  -h  d-h...; 


donc 


L’intégrale  générale  du  second  membre  est 


^ ^ ^ fi 

X "i  x-f-l  3 X 7. 


et  si  nous  la  prenons  entre  les  limites  o et  i , nous  pourrons  écrire 

d.lT[x\ 


dx 


\ x)  \2  x-hï  ) \3  ,r + 2 


d’où 


d\ir(x)  __  I 

dx‘  X’  (x  +i)»  ' (xH-2)-  “(,r-+-3)^ 


I I 


7 


( 5o  ) 


Considérons  actuellement  l’expression 


] =2“'»- 


En  développant  en  série  Z -4-  ^ ^ j » on  voit  que  u^.  sera  donné  pai  la 

formule 

U,n  - + 2 J [a;  + m ~ 2 (a:  + 3 (x  -4-  mp  J 

Prenons  deux  fois  la  dérivée  de  la  fonction  V,  il  vient 


+ /«  / 2 \ a:  H-  w 


-f- 


■ ...^1 

X + m + 1 / J 


et 


I t I ' r ^ I — 2 !• 

dx'^  ~~  2a\  X + m~^  x-\- ni+i'^  + (.r  + m4-i)Jl 

m—O  ^ 

Or  il  est  évident  que,  m recevant  toutes  les  valeurs  entières  de  o à l’infini, 
cette  dernière  expression  se  réduit  à 


d^Y 

dx^ 


• m=.^ 

~ X 2^  2 mY  1 


série  qui  sera  convergente  pour  toutes  les  valeurs  de  x plus  grandes 
que  O. 

Si  cette  série  est  convergente,  celle  qui  donne  et,  par  suite,  celle  qui 

donne  l’expression  de  V,  l’est  aussi. 

Cela  posé,  on  a 


I 

[x  -+■  Ifl) 


m=o 


dUV{x) 

d^^ 


Il  vient  alors 


ou 


donc 
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V I 1 ^ d^lr{x) 

dx'  X 2 x’  dx''  ’ 


d'‘lY{x) 

dx‘‘ 


I I rf’V 

■r  2x’  dx^  ■ 


d.lT(x 

dx 


) 1 ^ U 

- — Ix 1-  — -H  K, 

2 X dx 


et  en  intégrant  encore  une  fois, 


IT  [x)  = Yi  X X l X — X ~ - Ix  V: 

2 ’ 

de  là  résulte 


/r  ( 4- 1)  = A -h  -\-xlx  — x-\-~lx-\~y. 

2 


Il  reste  à déterminer  les  constantes  A et  B. 

La  formule  (a4)  (formule  de  Stirling)  est 

IT  ix  -Y-  \ ) = - l -in  X Ix  — X - Ix  • -5 

2 2 I .7.  X 

donc 


A — - lin  -Y-  Bj:  4-  V — 
2 


9B,  I 

1.2  X 


et  divisant  par  x, 


A — 4-/21 


B 


•r  L 1 .2  X J 


La  quantité  entre  crochets  se  réduit  à o pour  x = ^ 
nuis  pour  celte  valeur  de  x,  donc 


car  V et 


0B, 

I . 2 


B = O. 


1 

- sont 

X 


7- 


Cela  étant,  on  a aussi 


( 5a  ) 

A = - / 2 TT. 

2 

Donc  enfin  on  a la  série  convergente 

m=co 

(a5j  /r(^+ i)=^/27r+  [x-h^lx  — x+'^  u,n. 

• m=o 

Cette  formule  a été  donnée  pai’  M.  Liouville  dans  le  cours  qu’il  fait  au  Collège 
de  France.  J’en  dois  la  connaissance  à M.  Serret,  qui  a démontré  la  conver- 
gence de  la  série  ^ comme  je  l’ai  exposé  plus  haut. 


APPENDICE  RELATIF  AUX  INTÉGRALES  EULÉRIENTsTS  DE  PREMIÈRE 
ET  DE  SECONDE  ESPÈCE. 


22.  Je  terminerai  ce  travail  par  l’exposition  des  propriétés  des  intégrales 
eulériennes  dont  j’ai  eu  à me  servir. 

De  la  fonction  T {æ)  ou  intégrale  eulérienne  de  deuxieme  espèce. 

La  fonction  Ffjc)  est  caractérisée  par  l’intégrale  définie 

r(a:)  = C e~'^.cé^~'Kdcf., 

expression  dans  laquelle  x est  nécessairement  plus  grand  que  o.  Elle  n’anrait 
aucun  sens  pour  des  valeurs  nulles  ou  négatives  de  x. 

L’intégration  par  parties  donne 

r e~"-.a^^~'\da.  = [x — i)  j 

t'o  t/O 


(0 
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expression  dans  laquelle  x est  nécessairement  plus  grand  que  i . On  a donc 

r(a?)  = (x  — i).r(jc'  — i). 


Si  l’on  attribue  à x les  valeurs  entières 


4» • • • > 


il  vient 

r(a)  = r{,), 
F(3)  = ar(a), 
r(4)  = 3r(3), 


r^x)  = (x  — i)r(x-i). 

Donc 

r(x)  = 1.2.3. ..(æ*  -i)r(i); 

mais 

^00 

r (i)  = / e-'^da  = I, 


et,  par  conséquent, 


r {x)  — 1 .2.3. 4-. — i). 


25. 


CO 

f d a,  dans  laquelle  on  suppose  m > o. 

0 


Posant 


ma  — 


m 


et 


m 


d’où 


Il  vient 
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t/o 


et  remplaçant  m par  jc 


Jo 


Si  d’ailleurs  on  fait  n = i , 


r(«)  = r(i)  = i, 

et  l’on  arrive  à la  relation 


qui  est  évidente. 


2i,  Je  vais  maintenant  démontrer  la  formule  qui  m’a  servi  de  point  de 
départ  dans  le  n“  20.  On  a 


dx 

X 


J ^00 

I ^ .dx  . d a., 

0 

d’où,  intégrant  entre  les  limites  i etx, 


[a] 


Ix 


£ 


® ^ fiC  X 


•dcc. 


Dans  le  cas  où  x est  entier  on  peut  faire  successivement  x égal  à 
1 , 2,  3,...  {x  — i),  et  alors 


zr 


Xoo 

I (x  — i)  e~*  — [e~“  + 


d a 


..  ■+■  j 


et  comme 
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il  vient 

c”"  -f-  e~  2 » _ 

(26) 

-'-fi 

Jo  I t — e— * J a 


Cette  formule  (26)  n’est  ainsi  démontrée  que  pour  le  cas  de  x entier;  mais 
d après  le  point  de  vue  plus  général  sous  lequel  nous  avons  donné  le  développe- 
ment de  lT{x  +1),  il  est  indispensable  de  démontrer  qu’elle  est  vraie  sans 
restriction  pour  toutes  les  valeurs  positives  de  x. 

Voici  la  démonstration  adoptée  par  M.  Serret  dans  son  Algèbre  supérieure 
et  qui  ne  diffère  pas  essentiellement  de  celle  de  M.  Cauchy. 

Nous  avons 


r(x)=  Ç .da, 

t/o 

d où,  en  prenant  les  dérivées  par  rapport  à x, 

r'(x)=r  e-’^.a.^-' .la.da. 

Jo 

Mais  d’après  la  formule  (a), 

la=  f 

Jo 


d^-, 


donc 

rfxi  = 


J^co 

I e- a est  la  fonction  r(x), 
0 


(,S-hl)a  _ J , 


dj^ 

T' 


r(^) 

(P+0^’ 


Par  conséquent, 
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(■^) 

r(:r) 


(P+i)*J  P 


On  a d’ailleurs 


J dx  = e~é  — i)  — 


(P+Q-‘-(P+o-^ 

^(PH-0 


et  comme 

_ d.lr{x) 
r (j;)  dx  ’ 

on  peut  écrire 


/(p+i)  J P ■ 


Faisons,  dans  cette  expression,  x = 2 et  nous  aurons  en  vertu  de  la  relation 

r(2)=r(i), 

/r(2)  = zr(i)  = /.  1 = 0, 

d’où 


(P+i)-'-fP+Q-nrfp  _ . 

/(P  + i)  J P 


or 


donc 


et  aussi 


[(/3  + i)-‘  - (/3  + 0-]  = (/3  -F  (i-/3-0=-(/3  + 0-.  /3, 

Jo  r /(p+i)  J P - Jo  L P ^(p+«)J  ^ ’ 

r^{x-x)e-^  ,n_  p(^-l)(P+l)- 

Jo  . ^(P+0 


la  valeur  de  (x)  pourra  alors  se  mettre  sous  la  forme 

JY , s r”  I f/y  iWt  I ~i~ P)~‘ ~ + P)~'^1 

lT(x)  — j^  0(  +p)  P J- 


Posons 
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1 + ^ — d’où  /(i+|3)=a  et  </]3=e“ffa, 
les  limites  restent  les  mêmes,  et  il  vient 


^■^1  e*‘  d<x 

? 

a 


et  enfin 


/r(x)=  r(a:  — i)e~*  — ^ 1 

Jo  L J “ 

Ce  qui  est  la  formule  qu’il  s’agissait  d’établir  d’une  manière  générale. 


De  la  fonction  B q),  ou  intégrale  eulérienne  de  première  espece. 
Ses  relations  avec  la  fonction  r. 

25.  Cette  fonction  se  définit  par  l’égalité  suivante  : 


^{p.,q)=  f (i  — . dx., 

Jo 

où  X est  compris  entre  o et  1 , p e\.  q sont  positifs. 

Il  est  d’abord  facile  de  voir  que 


B [p.q)  = B(^  ./j), 

I 

c’est-à-dire  que  la  fonction  reste  la  même  quand  on  change  p q et  q en  p. 

En  effet,  dans  B(^,  p)  faisons  i — x = x' , d’où  x = i — x',  et  l’on  aura 
évidemment 


B (7»  P)=  f ■ (• 

, Jo 


xY~'  dx  — 


x'^P  . (i  — x')  dx'—  . . . = B (/?,  q). 


8 
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La  fonction  F peut  s’écrire 

T{p)=  Ç 

Jo 


et  en  posant 

il  vient 

(a) 

et  de  même 

{b) 


a = x^,  d’où  da='2xdx, 

r (p)  = a r . x^P~'  . dx, 
Jo 


F f g)  = 2 r e~^' . . dy. 

Reprenons  maintenant  l’expression 

B (p,  q)=  Ç ocP-'  (i  — x)’-'  . dx, 

Jo 

dans  laquelle,  avons-nous  dit,  x est  compris  entre  o et  i . Cela  étant,  soient 
x = singea,  i — x = cos“  w,  r/x  = 2 sin  w cos  wr/w. 

Pour  x=  O,  w=  O,  et  pour  x = i , w = donc 

7T  TT 

[d)  B(p,ç)=2  / sin*^“' w . cos^^"' uf/oj  = 2 Ç cos'*/’“'^  ca  . sin<^^~'' ca  r/w. 

^ c/o  t/Q 

Multipliant  les  égalités  {a)  et  (6)  membre  à membre,  il  vient 

F (p)  F (<y)  = 4 r • x^^“* . dx  Ç . dj 

t/O  t/0 

t/O  t/O 

et  si  nous  faisons 

(O 


( 5g  ) 

lunégrale  double  représentera  un  volunte  s’étendant  à 

l'innni  du  côte^  des  a.  et  des  j positifs,  compris  dans  le  trièdre  OXYZ  et  l.mité 
par  la  surface  donnée  par  l’équation  [c). 

Ma, s au  l,eu  d'estimer  ce  volume  en  employant  les  coordonnées  rectibenes 
on  peut  employer  les  coordonnées  polaires  dans  le  plan  des  XY  ^ ’ 

Posons  ^ 

Jr=pcosw,  j^=|3sinu. 

Le  volume  d’un  élément  prismatique  sera  dans  ce  système  représenté  par 

pzdp  do), 

r.ntégraje  relative  à p étant  prise  de  o à l’inBni,  et  l’intégivale  relative  à „ 
O à -•  On  a d’ailleurs 


x^p-'  =z  p^P-*  cos=^-‘  W, 

el  en  faisant  les  substitutions 
Ce  qui  peut  s’écrire  aussi 

7T 

= cosf*P-Ou.sin(*?-o«.  Ç" e~p^p‘^P^'^v-*  dp 


Or 


e-p^  pi g-,  dp  ^ ^ 

On  a d’ailleurs  par  la  formule  [d) 


^ O 


e-p\p^<P^'^-'Kd(p^)  = F( 


P 7)- 


d’où 


aT  cos^P-Vojsin^^-'.wi/w  = Br 

Jo 


P^i) 


r(p)-r(<7)  = B(p,  q).v{p  -t-  9j, 


8. 


ou  encore 


( 6o  ) 


(27) 


B{p,q)  = 


lS£bIlîl' 


De  là  on  tire  les  relations  suivantes: 

Si  q = i — p,  ce  qui  exige  qu’on  ait  /?  < i , on  aura 


2*^.  Si  P = 1 — P 


— 9 
2 


mais 


3°.  Donc  aussi  F = s/n. 


dcD  = n. 


27 . Reprenons  la  formule 

B{p,q)=  £ — 

dans  laquelle  nous  faisons 


\ — x= » d ou  dx  = , — ; — - 

1-t-r  ('H-J) 


Les  limites  de  l’intégrale  deviennent  o et  oo  ; soit  d’ailleurs  ^ = i — p,  on  aura 

B(p. 


Posons 


y — " î 

dy'—%  dz  et  jP~'  = z^" 


d’où 


il  viendra 


( ) 


B 


{p,  i-p)~  %n=  f £11 
J O >+2’ 


dz. 


Posons  aussi 


an/l  — 1 = 2 w,  d’où  p=. 


2 M + I 


2 n 


our  qu  une  va- 


I 

Ce  qui  est  permis,  car  on  peut  toujours  supposer  n assez  grand  p 

leur  / attribuée  à ^ diffère  de  d'une  quantité  aussi  petite  qu'on  voudra. 
D après  cela,  on  peut  écrire 

^{p,  p)  = 2n  f - 
do  * 

Mais  l’intégrale  étant  mise  sous  cette  forme,  on  voit  que  la  fonction 
change  pas  quand  on  change  z en  - z,  et  par  conséquent 

B{p,  i-p)  = r[p).r{i-p)  = 


dz. 


I-t-  2* 


ne 


2’' 


dz. 


Il  s'agit  de  trouver  cette  intégrale  définie,  et  pour  cela  considérons  la 


vante 


sui- 


dx. 


Les  racines  de  l’équation  binôme 


t H-  x^"—  O 


sont  de  la  forme 


cosp  ± ÿ'  — I sin  e, 

expression  dans  laquelle  u = «t  où  k ,*ut  prendre  toutes  les 

leurs  entières  de  o à « — i inclusivement. 

On  aura  donc 


va- 


f ~.rdx=  V ' 

— COS(t)’+ 

A=o 


( 62  ) 


L’intégrale  générale  / , 

^ ^ J \x  — cosv 


B 


l [(j:  — cosp)*  + sin^t^] 


)^-f-  sin’»» 

A cosp  + B 


dx  est  à la  constante  près 

^ [ X — cos  V \ 

arc  tang  - — ^ ) 

^ \ sin  I'  ) 


sine 


donc 


A‘=n — I 


^ A ^ r(x  — cosp)’-H  sin’eT 
2 L(x  4- cose)^ -J- sin’ e J 


k—o 

k=zn — I 


A cos  e + B 


k=o 


t(  X — cose\  /Æ-4-cose 

“'■‘=  ^iSr-  ) + ‘fc  "g  ( -lis;;- 


et  pour  les  limites  de  — qo  et  4-  co  , il  vient 


A=rn--i 


r“_iüL^2  = ;r  y 

/ i4-  z’"  ^ 

•^—00  .T, ^ 


A cos  e -t-  B 


La  décomposition  des  fonctions  rationnelles  donne 


d’où 


cos  ( 2 »2  -f-  1 ) e cos  2 mv 

A = — — — 5 ü — ' ■ ■ ; 


A cos  e -(-  B 1 ■ / \ 

= - sin  (am  4-  i)  v. 


sm  V n 


Mais  on  a écrit  plus  haut 


‘î.m  \ = inp  et  \>  = 


[ik  i)7t 


2 n 


donc 


et,  par  suite, 


A cos  e 4-  B i . , , , 

= - sin  (a  A‘  4-  I j ^7:, 


sin  V n 


k = {n—i) 


= ï 2 sin(-2i  + .: 


pn. 


A = 0 


En  vertu  de  la  formule 


(63) 


siiirt  + sin(rt  -h  b)  + sin(a  -h  2b).  . . + sin[rt  + («-  i)^] 

2(1  — cos  F)  ' ’ 


nous  aurons 


y ^\n{<ik i) pu  = 

, I cos  2 B 7T 

* = O ' 


sin  /3  TT  I 

sin’/>77  sin/?7r’ 

donc 


sin/j7r[i  — cos(2  w + i)it] 
I — cos  2 TT 


r(p).T(x-p)  = 


r”  z'” 


dz  = 


sin  pn 


1 


ou,  ce  qui  est  la  même  chose, 


28.  Arrivons  enfin  à la  formule  {a),  u°  20,  qui  nous  a servi  dans  la  for- 
mule de  Stirling. 

Prenant  les  logarithmes  dans  la  formule  (28),  il  vient 
' /r(jr)H-/I(f — x)  = lu  — l sinujc, 

et  aussi 


J'  djc -h  J'  ir  (î x)  .dx  = lu  — Ç l.sinux.dx. 

” ° Jo 


Op  voit  sans  peine  que  les  différentes  valeurs  par  lesquelles  passe  la  fonc- 
tion r(x)  pour  les  valeurs  de  x comprises  entre  o et  i sont  égales  à celles 
par  lesquelles  passe  la  fonction  r(i  — x). 

Donc 


2 f r(.r)  = lu  — / sin  7: JC, 


( 64  ) 


et 


Jf  IV [x)  dx  :=-ln  — - f l.sinnx. 

0 


Soit  maintenant 


I a o> 

nx  = co^  clx  = — 9 


les  limites  de  l’intégrale  se  changent  en  o et  -k. 
Et  il  vient 


Jr  lT(^x)dx  — -ln j Isittudu. 

0 


Mais  de  « à 7T  les  sinus  repassent  par  les  mêmes  valeurs  ; donc 


f Z.sin  w.f/w  = 2 r / sin  co.<^w  = 2 r l.coscô.dc^. 

»/0  c/0 


La  relation 


donne 


d’où  l’on  tire 


1 I 

sin  CO  = 2 sin  - co.cos-  co, 

2 2 ' 


/ sin  CO  = / 2 4-  / sin  - CO  4-  Z cos  - co  ; 


1 


. I dto  n J ^ ' 

Z.sin co.cZco  = 7r/2  4- 2 J 2 ^ J 

7T  ^ 

= r:/2  + 2 jf4sinM.rfo)  + 2 r I cos  co.cZco 
= 7.  J'  l sin  co.i/co; 


par  conséquent 


par  conséquent 

( 65  ) 

1 l sin  «f/w  = — 7r/2, 

ou  bien 

* \ 

«/O  ^ 2 

Remarquons  aussi  que 


il  vient  alors 

II 

1 

f + l.x]c/a:  = - l2n  — 

«/o  2 . 


D’ailleurs 

et  puisque 

“t-  ix  — /xr(x), 

on  peut  écrire 

xr(x)  = r(x-i-j}, 

f'  l-Tfx  ~h  i).dx  = - /an  — i, 

Jo  2 > 

qui  est  la  formule  qu’il  s’agissait  de  démontrer;  elle  est  due  à M.  Raabe. 


9 


( 66  ) 


PROGRAMME  DES  QUESTIONS  D’ASTRONOMIE. 

1 . Notions  historiques  sur  la  découverte  de  la  gravitation  universelle. 

2.  Recherche  de  la  cause  qui  fait  décrire  aux  planètes  et  à leurs  satellites 
des  orbites  elliptiques. 

Lois  de  Répler;  conséquences  qui  en  dérivent. 

Lois  de  la  gravitation  pour  les  planètes  et  les  satellites;  cas  particulier  de 
la  lune,  satellite  unique  de  la  terre. 

5.  Les  corps  célestes  agissent  les  uns  sur  les  autres,  en  vertu  de  leur  forme 
sensiblement  sphérique,  comme  s’ils  étaient  réduits  à des  points  mathéma- 
tiques; attraction  de  deux  sphères  ou  de  deux  couches  sphériques  homo- 
gènes. 

i.  Masse  des  planètes  comparée  à celle  du  soleil;  méthode  spéciale  pour 
la  terre. 

r>.  Mouvement  d’un  point  attiré  vers  un  centre  fixe  en  raison  inverse  du 
carré  de  la  distance. 

(>.  Formules  du  mouvement  elliptique. 

7 Méthode  de  Lagrange  pour  le  développement  de  certaines  fonctions 
implicites. 

8.  Solution  du  problème  de  Képler. 

Fu  et  approuvé^ 

Strasbourg,  le  ag  août  1857. 

Le  Doyen, 

L.  DAUBRÉE. 

Permis  cl  imprimer, 
Strasbourg,  le  3o  août  1857. 

Le  Recteur, 

DELCASSO. 
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